Sharing Social Network Data: Differentially Private Estimation of Exponential-Family Random Graph Models

Motivated by a real-life problem of sharing social network data that contain sensitive personal information, we propose a novel approach to release and analyze synthetic graphs in order to protect privacy of individual relationships captured by the social network while maintaining the validity of statistical results. A case study using a version of the Enron e-mail corpus dataset demonstrates the application and usefulness of the proposed techniques in solving the challenging problem of maintaining privacy \emph{and} supporting open access to network data to ensure reproducibility of existing studies and discovering new scientific insights that can be obtained by analyzing such data. We use a simple yet effective randomized response mechanism to generate synthetic networks under $\epsilon$-edge differential privacy, and then use likelihood based inference for missing data and Markov chain Monte Carlo techniques to fit exponential-family random graph models to the generated synthetic networks.Comment: Updated, 39 page

Differentially Private Model Selection with Penalized and Constrained Likelihood

In statistical disclosure control, the goal of data analysis is twofold: The released information must provide accurate and useful statistics about the underlying population of interest, while minimizing the potential for an individual record to be identified. In recent years, the notion of differential privacy has received much attention in theoretical computer science, machine learning, and statistics. It provides a rigorous and strong notion of protection for individuals' sensitive information. A fundamental question is how to incorporate differential privacy into traditional statistical inference procedures. In this paper we study model selection in multivariate linear regression under the constraint of differential privacy. We show that model selection procedures based on penalized least squares or likelihood can be made differentially private by a combination of regularization and randomization, and propose two algorithms to do so. We show that our private procedures are consistent under essentially the same conditions as the corresponding non-private procedures. We also find that under differential privacy, the procedure becomes more sensitive to the tuning parameters. We illustrate and evaluate our method using simulation studies and two real data examples

Statistical Issues in Assessing Hospital Performance

From the Preface: The Centers for Medicare and Medicaid Services (CMS), through a subcontract with Yale New Haven Health Services Corporation, Center for Outcomes Research and Evaluation (YNHHSC/CORE), is supporting a committee appointed by the Committee of Presidents of Statistical Societies (COPSS) to address statistical issues identified by the CMS and stakeholders about CMS’s approach to modeling hospital quality based on outcomes. In the spring of 2011, with the direct support of YNHHSC/ CORE, COPSS formed a committee comprised of one member from each of its constituent societies, a chair, and a staff member from the American Statistical Association, and held a preliminary meeting in April. In June, YNHHSC/CORE executed a subcontract with COPSS under its CMS contract to support the development of a White Paper on statistical modeling. Specifically, YNHHSC/CORE contracted with COPSS to “provide guidance on statistical approaches . . .when estimating performance metrics,” and “consider and discuss concerns commonly raised by stakeholders (hospitals, consumer, and insurers) about the use of “hierarchical generalized linear models in profiling hospital quality. The committee convened in June and August of 2011, and exchanged a wide variety of materials. To ensure the committee’s independence, YNHHSC/CORE did not comment on the white paper findings, and CMS pre-cleared COPSS’ publication of an academic manuscript based on the White Paper

Statistical Inference in a Directed Network Model with Covariates

Networks are often characterized by node heterogeneity for which nodes exhibit different degrees of interaction and link homophily for which nodes sharing common features tend to associate with each other. In this paper, we propose a new directed network model to capture the former via node-specific parametrization and the latter by incorporating covariates. In particular, this model quantifies the extent of heterogeneity in terms of outgoingness and incomingness of each node by different parameters, thus allowing the number of heterogeneity parameters to be twice the number of nodes. We study the maximum likelihood estimation of the model and establish the uniform consistency and asymptotic normality of the resulting estimators. Numerical studies demonstrate our theoretical findings and a data analysis confirms the usefulness of our model.Comment: 29 pages. minor revisio

Test of candidate light distributors for the muon (g−-2) laser calibration system

The new muon (g-2) experiment E989 at Fermilab will be equipped with a laser calibration system for all the 1296 channels of the calorimeters. An integrating sphere and an alternative system based on an engineered diffuser have been considered as possible light distributors for the experiment. We present here a detailed comparison of the two based on temporal response, spatial uniformity, transmittance and time stability.Comment: accepted to Nucl.Instrum.Meth.

Gravitational and electromagnetic fields of a charged tachyon

An axially symmetric exact solution of the Einstein-Maxwell equations is obtained and is interpreted to give the gravitational and electromagnetic fields of a charged tachyon. Switching off the charge parameter yields the solution for the uncharged tachyon which was earlier obtained by Vaidya. The null surfaces for the charged tachyon are discussed.Comment: 8 pages, LaTex, To appear in Pramana- J. Physic

R.A.Fisher, design theory, and the Indian connection

Design Theory, a branch of mathematics, was born out of the experimental statistics research of the population geneticist R. A. Fisher and of Indian mathematical statisticians in the 1930s. The field combines elements of combinatorics, finite projective geometries, Latin squares, and a variety of further mathematical structures, brought together in surprising ways. This essay will present these structures and ideas as well as how the field came together, in itself an interesting story.Comment: 11 pages, 3 figure

Dopamine D1 vs D5 receptor-dependent induction of seizures in relation to DARPP-32, ERK1/2 and GluR1-AMPA signalling.

Recent reports have shown that the selective dopamine D(1)-like agonist SKF 83822 [which stimulates adenylate cyclase, but not phospholipase C] induces prominent behavioral seizures in mice, whereas its benzazepine congener SKF 83959 [which stimulates phospholipase C, but not adenylate cyclase] does not. To investigate the relative involvement of D(1) vs D(5) receptors in mediating seizures, ethological behavioral topography and cortical EEGs were recorded in D(1), D(5) and DARPP-32 knockout mice in response to a convulsant dose of SKF 83822. SKF 83822-induced behavioral and EEG seizures were gene dose-dependently abolished in D(1) knockouts. In both heterozygous and homozygous D(5) knockouts, the latency to first seizure was significantly increased and total EEG seizures were reduced relative to wild-types. The majority (60%) of homozygous DARPP-32 knockouts did not have seizures; of those having seizures (40%), the latency to first seizure was significantly increased and the number of high amplitude, high frequency polyspike EEG events was reduced. In addition, immunoblotting was performed to investigate downstream intracellular signalling mechanisms at D(1)-like receptors following challenge with SKF 83822 and SKF 83959. In wild-types administered SKF 83822, levels of ERK1/2 and GluR1 AMPA receptor phosphorylation increased two-fold in both the striatum and hippocampus; in striatal slices DARPP-32 phosphorylation at Thr34 increased five-fold relative to vehicle-treated controls. These findings indicate that D(1), and to a lesser extent D(5), receptor coupling to DARPP-32, ERK1/2 and glutamatergic signalling is involved in mediating the convulsant effects of SKF 83822

Design and Performance of SiPM-Based Readout of PbF\u3csub\u3e2\u3c/sub\u3e Crystals for High-Rate, Precision Timing Applications

We have developed a custom amplifier board coupled to a large-format 16-channel Hamamatsu silicon photomultiplier device for use as the light sensor for the electromagnetic calorimeters in the Muon g - 2 experiment at Fermilab. The calorimeter absorber is an array of lead-fluoride crystals, which produces short-duration Cherenkov light. The detector sits in the high magnetic field of the muon storage ring. The SiPMs selected, and their accompanying custom electronics, must preserve the short pulse shape, have high quantum efficiency, be non-magnetic, exhibit gain stability under varying rate conditions, and cover a fairly large fraction of the crystal exit surface area. We describe an optimized design that employs the new-generation of thru-silicon via devices. The performance is documented in a series of bench and beam tests